debug=1;
c=1480;
vel=Doppler2D(ans(:,500:800)',40e6,5e6,0.6,2e3,1540,debug);

%%
iq=pwdoppler;%(:,1:14:end);
%iqz=[iq(1:15) zeros(1,43)]
%iq2=[iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iq;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz;iqz];
vel=Doppler2D(iq(:,39:58)',40e6,5e6,0.6,2e3,1540,debug);
%%
P=pwdoppler';
%%
p=P(800:-1:500,:);
dec=ceil(fs/(B*fc));
p=p(1:dec:end,:);

fs=40e6;
fc=5e6;
debug=1;
c=1540;
prf=20000;
B=0.6;
vel=Doppler2D(p,fs,fc,B,prf,c,debug);
%%
vel=tand(22.8)*(22/38)*(1540/2)*(dec/fs)*prf

%tand(26.8)

%tand(14.7)